# • 



Mass Spectrometer ; 1 Computing Device 



,1 15 ^ 



FIG. 1 



# 



L RECEIVE DATA 



DENOISE 



Wavelet Transformation 



Noise Profile 



Apply Noise Profile to each Stage 



f>0 



BASELINE 
CORRECTION 



Generate Sparse Data Set 



*~\ i ■ n /~v ■ i 

ornii oiyiidi 



Generate Peak-Free Signal 

Generate Baseline 

Subtract Baseline From Signal 



COMPRESSION 



PEAK FITTING 



HFNOTYPF 



Mass Shift 

Find width of Strongest Peak 
Delete an Area Around Each Peak 
Generate a Residual Baseline 
Subtract Residual Baseline From Signal 
Determine Height of Each Peak 

Determine Peak Probabilities 



FIG. 2 



50*1- 




4500- 




4000- 


3500- 




^3000- 




| ?M~ 




~ 2000- 




1500 -j 


looo-i 

1 


500- 




111 


1 

-3S2-, 



fill 



lit . 



H« ...... iK H'&J-Ji 



^-4u • - . y, I . ; 

' *** . ™~ -1"*^. \_ lie ^i, ■ ■ - ■ - -> 



7000 8000 SOGO 10099 10753. 



Mm* 



FIG.3 



Mass 
Spec. 
Data 



J 

no 



Stage 0 
Hi 



I Stage 0 j 
* Lo f 



Stage 1 
Hi 



Stage 1 
Lo 



8-3 



HStage 2 
* 1 Hi 

I Stage 2 
*i Lo 



8 1, 



0^ 



12- 

IT 



0 



Stage n : 

__Hi 

Stage n 
Lo 



FIG. 4 




FIG. 6 



1868- 

-11 48-? 
0 



1806- 

■1790-, 
0 



1st Stjg* 



High Pass \ \L J ' t 



2900 5000 7500 1020a 
3p0 3tJe* Hl«h fm 



2000 256?: 



Cc#npf**aon iggO- 

1.00- .. -I945-? 

0 

aao- " 



1805- 
4152 



1000 

ijtfc Slag* Hijh P*» - 

fl ( »t ^ m *ti> vi i i i»M *M ii> H v ^ 0.20- 

0 200^f 400 B5^ ft()o: 




2000 



4000 5111 



1382- 
•193S- 



4th Stij« High Pjb 



21231- 



0 500 1000 1295 

01h Stj«« U»?w 




FIG. 7 



Stage 0 - Hi 




Threshold 0=4XNoiseProfile 




Stage 1 - Hi 



Threshold 1=2XNoiseProfile 




Stage 2 - Hi 



Threshold 2=1XNoiseProfile 



PS 




Stage n - Hi 



Threshold n=(1/2 n - 2 )XNoiseProfile 



( I 
r'S 



Stage n - Lo 



FIG. 8 



1rf SUg« High P« 



148- 

8 2500 MOO 7500 181 », 



2n4 Stif • Miflh Pi» 



548- 



3rd Stif • Ifcgfc P*s 



-413 , 
0 

2829-. 



*+ f * 4 



■I : 



1000 2000 2567. 



344$-, 
0 



206 



400 



653! 



1 00- 
080^ 

040-' 

020^ 

ooo : 



459-,' 
0 

869- 



2000 



4000 5111 



1024-, 
21488 



4th SU«« HijS Pj«* 




500 1000 123S 

*44h Staff LwiPjai 




FIG 9 



Sj g na l (t)= (Start Oft) + Start 1 (t) + Start 2 m... + Start 23 ft)) 

FIG. 10 - SHIFT SIGNAL TO ACCOUNT FOR 
VARIATIONS DUE TO STARTING POINT 



Stage 0 

Stage 1 — LpJ LL 

Stage 2 Lj_J — 

Stage n-hi I f 1 I 



FIG 11 




309 



1 

1*» 



4 



I > 



1003* 
503* 



1 1 .. 
: 1 . , r 



300Q 4000 MOO (000 7000 1000 



9000 100(0 1973) 



FIG 13 - TAKE A MOVING AVERAGE, REMOVE SECTIONS EXCEDING A 

THRESHOLD 




i — i /—\ a a r— i 



nu. \H - rilNU IVI l IM I IVI/A IIN KtlVIAIININlj £>IUNALb ANU UUNIMbU I IU 
FORM A PEAK FREE SIGNAL 



DO A WAVELET 
TRANSFORMATION ON 
THE PEAK FREE SIGNAL 



111 



DELETE ALL STAGES, ' 
LEAVING ON THE N-LO !— | 
STAGE 



RECONSTRUCT A 
SIGNAL FROM N-LO 



GENERATE BASELINE 
CORRECTION SIGNAL 



FIG. 15 - GENERATE BASELINE CORRECTION 




FIG 16 




eooc 



7000 



9000 



SOB 



10000 



11322 



FIG. 17 



Mais 



NON-0 

COEFFICIENTS 



VALUE 



100 
150 
500 
10,050 
10,075 
11,125 
12,100 
13,250 



IS 



25 
220 

.1 

800 
890 
910 
1000 (MAX) 
940 



ft 



V IB* 

I NTERMEDIAT E 

100.025 
150.220 
500.0001 
10,050.8 
10,075.89 
11,125.91 
12,100.99999 
13,250.94 




RELATIVE 

100.025 
50.220 
350.0001 

9550.8 
25.89 
150.91 
975.99999 

1150.94 



FIG. 18 



'if 



RECEIVE DATA FILE 



10* 




ARRANGE AN X-Y ARRAY 
OF COEFFICIENTS 



16 3 



DETERMINE STARTING 
POINT AND MAXIMUM 
DATA VALUE 



GENERATE 
INTERMEDIATE REAL 
NUMBER 



CALCULATE 
COMPRESSED VALUE 



Whole Portion is X axis number 

Decimal Portion is fractional Data 
Value 



FIG. 19 





FIG. 21 



A . A Z 



FIG. 25 




FIG 26 




FIG. 28 



32 S 




<ftfL 



33/ 



PROBABILITY OF GG EXISTING: 

P(GG) = P(G) *P(1-C) 

=90% * (100% -20%) 
=90% * 80% 
=72% 



PROBABILITY OF CC EXISTING: 

P(CC) = P(C)*P(1-G) 

=20% * (100% -90%) 
=20% * 10% 



= 2% 



333 



PROBABILITY OF GC EXISTING: 

P(GC) = P(G) *P(C) 
=90% *20% 
= 18% 



FIG. 31 



3 



3^' 



J- 



FIND H OF EACH PEAK ^ 

\ V 

EXTRAPOLATE NOISE PROFILE 

FIND NOISE OF EACH PEAK 

CALCULATE A SIGNAL TO NOISE RATIO *7 ' 

FIND A RESIDUAL ERROR — 1 

CALCULATE AN ADJUSTED SIGNAL TO '1 Y 
NOISE RATIO 

DEVELOP A PROBABILITY PROFILE ^ 1 



DETERMINE PEAK PROBABILITIES 



I 

\K 

1 

AL 

I 



ADJUST PEAK PROBABILITY BY ALLELIC 
PENALTY 



CALCULATE GENOTYPE PROBABILITIES 



DETERMINE AN ALLELIC PENALTY 
1TY 

. PE PROBABILITIES ^ ^ ^ 

\ 

TEST RATIO OF GENOTYPE ~ ^ ^ 

PROBABILITIES 

.{ . ^6 

FIG 33 



